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Abstract 

Solid-solution strengthening results from solutes impeding the glide of dislocations. Existing theories of strength rely on solute- 
dislocation interactions, but do not consider dislocation core structures, which need an accurate treatment of chemical bonding. 
Here, we focus on strengthening of Mg, the lightest of all structural metals and a promising replacement for heavier steel and 
aluminum alloys. Elasticity theory, which is commonly used to predict the requisite solute-dislocation interaction energetics, is 
replaced with quantum-mechanical first-principles calculations to construct a predictive mesoscale model for solute strengthening 
of Mg. Results for 29 different solutes are displayed in a "strengthening design map" as a function of solute misfits that quantify 
volumetric strain and slip effects. Our strengthening model is validated with available experimental data for several solutes, includ- 
. ing Al and Zn, the two most common solutes in Mg. These new results highlight the ability of quantum-mechanical first-principles 
calculations to predict complex material properties such as strength. 
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1. Introduction 

Inherent limitations of strength and formability, which are 
related to microscopic-scale deformation mechanisms, are sig- 
nificant obstacles to widespread adoption of wrought magne- 
sium in transportation industries. Magnesium's poor workabil- 
ity at room temperature comes from the anisotropic response 
of its hexagonal closed-packed (HCP) crystal structure, and the 
strength of conventional Mg alloys is lower than that of most 
aluminum alloys without added processing steps (e.g. grain 
refinement). While the polycrystalline ductility of face- and 
body-centered cubic metals results from multiple symmetry- 
related slip systems, the basal and prismatic planes in Mg are 
perpendicular in its HCP lattice, unrelated by symmetry, and 
must both be active to achieve appreciable ductility. The room 
temperature stress required to plastically deform Mg along its 
(easy) basal slip plane is two orders-of-magnitude lower that 
the (hard) prismatic plane. The five independent slip systems 
required by the von Mises criterion|[T] for sufficient ductility 
are simultaneously activated only at temperatures near 300°C. 
Insights required to overcome these design challenges can come 
from new predictive capabilities. To that end, we develop a new 
accurate first-principles strengthening model that predicts Mg 
basal strengthening as a function of substitutional solute chem- 
istry via computation of the fundamental solute-dislocation in- 
teraction. 

We compute the interaction energy of solutes with screw 
and edge basal dislocations — finding surprisingly similar mag- 
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nitude interactions of solutes with both dislocations types — and 
use this information to build a first-principles "design map" for 
the strengthening of solutes in Mg in a computationally efficient 
manner. Using chemically accurate predictions of the atomic- 
scale geometry, we resolve the volumetric expansion and com- 
pression and the local slip in both the far-field and into the dis- 
location cores, and determine the solute-dislocation interaction 
energy for all substitution sites. This combines accurate disloc- 
ation core geometries at the atomic scale — available from state- 
of-the-art first-principles calculations|2][3]|4][5]|6) with flexible 
boundary condition methods[7, 8 1, proven successful for Mo|9| 
and Al|10| — and computation of solute "misfits." The misfit 
of a solute quantifies relative changes in the lattice from dilute 
substitution of a solute, and consists of two important compo- 
nents: a "size" misfit for the change in the local volume, and 
a "chemical" misfit for the change in energy to slip the crystal 
in the basal plane. The misfits provide the basis for approx- 
imating the solute-dislocation interaction energy; we test this 
approximation against direct substitution of Al in Mg disloc- 
ations before applying it for other solutes. The solute/screw 
dislocation interaction energy, normally assumed to be negligi- 
ble in the far-field from elasticity treatments lfTTI fT2l but which 
has recently been quantified with atomic-scale studies(9] [T3l , 
is found to be nearly identical in magnitude to the edge dis- 
location interaction energy due to the screw core geometry and 
contributes to the prediction of strength. Finally, we use our 
geometrically-informed calculation of solute-dislocation inter- 
actions to predict the dilute-concentration strengthening effect 
of 29 different solutes and create a simple map of strengthen- 
ing potencies that suggests new Mg alloy designs. Our entirely 
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first-principles approach is validated with available experimen- 
tal strengthening data. 

2. Computational methods 
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Figure 1 : First-principles phonon spectra (top) and (0001) generalized stacking- 
fault energy (bottom) for Mg. These describe response of the lattice to small 
displacements, elastic deformation, and slip in the basal plane. Density- 
functional theory is able to accurately reproduce the vibrational spectra from 
experimental measurements: Neutron scattering data from |14|, and a Born- 
von Karman fit to the data from 1 15]. The equilibrium lattice constants from 
the ultrasoft pseudopotential are a = 3.19A and c = 5.18A; the elastic con- 
stants, which determine the slopes of the spectra for long wavelengths are 
Cn = 60GPa, C33 = 61GPa, Cn = 21GPa, C13 = 20GPa, and C44 = 18GPa. 
For the generalized-stacking fault surface, a single plane of Mg is displaced in 
the basal plane by a linear combination of j[l 120] and a[ll00]; the defected 
geometry is allowed to relax in the [0001] directions, and the energy per area 
for the defect is the generalized stacking-fault energy. The highlighted point of 
|[1010] is a metastable configuration known as the intrinsic 12 stacking fault. 
The 712 energy agrees well with other density-functional theory calculations of 
the same| 16 1. The 12 stacking fault geometry is the basis for computing chem- 
ical misfits for solutes. 

Efficient computation of interactions between solutes and 
dislocations requires simulation cells and flexible boundary con- 
dition approaches that isolate individual defects. Calculations 
are performed with VASP| 2 3|, a plane-wave density -functional 
theory code. Magnesium and all solutes are treated with Van- 
derbilt ultrasoft pseudopotentials|5, 6], and the Perdew-Wang 



91 GGA exchange-correlation potential|4|. The ultrasoft Mg 
pseudopotential ([Ne]3i 2 ) accurately reproduces experimental 
lattice! 17| (l ess than 0.9% error) and elastic] 18 1 (less than 5% 
error) constants, and phonon frequencies! 14, 15 1 (less than 3% 
error) of bulk Mg; c.f. Figure [T] The stacking-fault surface 
in Figure [T] while not available experimentally, compares well 
with other density-functional theory calculations [ 1 6 1 . We chose 
a planewave cutoff of 138eV, fc-point meshes (see below for 
specific values tied to each geometry) and Methfessel-Paxton 
smearing of 0.5eV to give an energy accuracy of 5meV for 
bulk Mg. For calculations involving solutes, the cutoff energy 
was increased as necessary to accommodate harder pseudopo- 
tentials for solutes; c.f. Table [2] for all cutoff energies used for 
the corresponding pseudopotentials (in general, we selected a 
cutoff energy of 1.3 times the suggested cutoff for the poten- 
tial). The geometries for misfits, and dislocation calculations 
and coupling with lattice Green function flexible boundary con- 
dition methods are given below. 

2.1. Size misfit 

For the size misfit, we substituted single solutes into a 2 x 
2x2 Mg supercell (£-point mesh of 16 x 16 x 10) at five differ- 
ent volumes based on the equilibrium Mg volume Vq: I.I6V0, 
1.05V , l.OOVo, 0.95V , 0.86V . The atomic positions in each 
supercell were relaxed until all forces were less than 5meV/A. 
The size misfit is the logarithmic derivative of the Burgers vec- 
tor b — I [21 10] with solute concentration c s in the dilute limit; 
hence, 
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The change in energy for a solute with volume is given by the 
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where B is the bulk modulus. We use the slope £^ olllte to com- 
pute the change in binding energy for a solute due to local 
volumetric strain. Comparison with calculations using a larger 
3x3x3 supercell affected the size misfit by < 10%. 

2.2. Chemical misfit 

For the chemical misfit, we substituted single solutes into a 
2x2x9 Mg supercell (£-point mesh of 16x16x1) correspond- 
ing to 18 (0002) planes. The supercell vector along the c[0001] 
direction has an extra |[0110] component, so that the supercell 
represents a periodic repetition of stable stacking faults sepa- 
rated by a distance of 9c. The solute is substituted into a site 
in the stacking fault, and the atomic positions were relaxed 
until forces were less than 5meV/A. Because the 12 intrinsic 
stacking-fault configuration is (meta)stable (c.f. Figure fTJ, we 
do not impose any constraints on the relaxation. The chemical 



misfit is the logarithmic derivative of the 12 intrinsic stacking- 
fault energy jn with solute concentration c s in the dilute limit; 
hence, 
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The chemical misfit is calculated from our supercell energies as 
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where the "displaced" and "undisplaced" geometries correspond 
to the layered structure with and without an 12 intrinsic stacking 
fault, and V3a 2 /2 is the basal plane area. 

2.3. Dislocation geometries 
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Figure 2: Full Mg basal screw (top) and edge (bottom) dislocation core equi- 
librium geometries, separated into regions I (blue), II (red), and III (gray) and 
identified solute sites (green). The periodic simulation boxes are 69. 14x54. 42x 
3.19 A 3 for screw and 84.92x72.57x5.53 A 3 for edge. The initial screw geom- 
etry has displacements coming out of the page, while the edge displacements 
are in the plane of the page; relaxation produces additional displacements in the 
core of the partial dislocations. 

Flexible boundary condition methods lfT9H 71[8l relax the pure 
Mg basal dislocation geometries. The starting geometry comes 



from the anisotropic elasticity solution for a dislocation using 
the elastic and lattice constants from our first principles calcu- 
lations COl |2D . The elasticity solution determines the displace- 
ment field as a continuum function; we start with a bulk HCP 
lattice periodically repeated along the dislocation line direction: 
| [2110] for the screw (fc-point mesh of 1 x 1 x 16) and o[0110] 
for the edge (fe-point mesh of 1 x 1 x 12). The displacement 
field is centered between two (0002) atomic planes, and applied 
to every atom. A finite sized computational cell is produced by 
simulating only atoms within 10 lattice planes from the esti- 
mated positions of the partial cores; c.f. Figure[2] The distance 
from the core is determined by summing the distance to return 
to small relative displacements (the size of region I), the dis- 
tance for the lattice Green function to match the elastic Green 
function| 8 1 (the size of region II), and the distance from an un- 
relaxed free surface to produce zero forces (the size of region 
III). Together, this ensures that region II both does not have spu- 
rious forces due to the vacuum outside and is sufficiently sepa- 
rated from the defect core to follow the bulk harmonic response. 
In the case of the screw dislocation, the lack of a long-range vol- 
umetric strain field allows the use of a periodic simulation cell 
with "domain boundaries" at the cell boundaries. For the edge 
dislocation, vacuum is needed to isolate periodic images of the 
dislocations. 

2.4. Flexible boundary conditions 

Lattice Green function-based flexible boundary conditions 
isolate the dislocation at the center of the simulation cell from 
the domain boundaries or vacuum while displacing surround- 
ing atoms as if the dislocation were embedded in infinite bulk 
responding harmonically |[T9l 171 . The method relaxes each atom 
based on its location within one of three different regions deter- 
mined by the distance from the dislocation core. Overall, the 
screw dislocation consists of 525 atoms and the edge disloca- 
tion consists 806 atoms. Region I atoms near the core (54 for 
the screw, 130 for the edge) start the simulation with non-zero 
forces, and are relaxed using a conjugate-gradient method. As 
the relaxation commences, the neighboring atoms in region II 
(164 for the screw, 247 for the edge) start with zero forces but 
the forces increase as atoms in region I are displaced. The re- 
gion II atoms can be treated as if they were coupled to infinite 
harmonic bulk; hence, their forces can be relaxed by applying 
a displacement to each atom R given by the lattice Green func- 
tion (LGF): u(R) = £# G(R - R')f(R'), where R' only varies 
over atoms in region II, / are their forces, and G(R - R') is 
the (tensor) LGF. This displaces all atoms — including region I 
atoms, and the outer region III atoms (307 for the screw, and 
429 for the edge). The atoms in the outer region have non-zero 
forces due to the domain boundary or vacuum; but the lattice 
Green function displaces them as if they were part of an infi- 
nite bulk lattice. The relaxation cycle continues between region 
I (conjugate gradient) and II (lattice Green function) until the 
forces are less than 5meV/A in both regions. The final result 
of the relaxation is the stress-free dislocation core equilibrium 
geometry. 



2.5. Direct solute-dislocation interaction energy calculation 

We compute the interaction of Al with Mg dislocations by 
direct substitution of Al for Mg at different sites in the disloca- 
tion cores; c.f. Figured For each substitution, the entire region 
I was relaxed until the forces were less than 5meV/A. This de- 
fines the relative energies for Al in each site (18 for the screw, 
30 for the edge); to define the energy zero for Al — the refer- 
ence of Al substituted into bulk Mg with no strain field — we 
reference the average energy of Al above and below the stack- 
ing fault region. The periodic repetition of the solute along the 
dislocation line introduces a finite-size error in the calculated 
solute-dislocation interaction. For the screw dislocation geom- 
etry, there is one Al atom every 3.19A along the dislocation 
line, and for the edge dislocation, every 5.53A. Using a screw 
dislocation geometry with double the periodicity (6.38A), the 
calculated Al solute/screw dislocation interaction energy was 
8.4meV higher than that extracted from the original geometry. 
We expect this to be an upper limit on the finite-size error as 
it is (a) for the shortest Al-Al repeat distance, and (b) for the 
largest change in local geometry. 

2.6. Misfit approximation of solute-dislocation interaction en- 
ergy 

To compute the interaction of any solute with the Mg disloc- 
ations, we analyzed each dislocation geometry in terms of local 
volumetric strain and slip. The local volumetric strain at each 
atomic site in the final relaxed dislocation geometry is defined 
from the nearest-neighbor positions as 
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where j? are the vectors to nearest neighbors for a site, and 
x are the corresponding nearest-neighbor vectors in the HCP 
lattice 1 22 1 . The slip interaction energy (E s n p ) is calculated at 
each atomic site as 
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where d are the vectors to the nearest neighbors in adjacent 
basal planes, y(oooi) is the generalized basal stacking-fault en- 
ergy for displacement d, and the factor of i is from assigning 
half the bond energy for the 3 out-of-plane neighbors. The in- 
teraction energy for a solute is then a sum of two contributions: 
the size interaction (given by the change in solute energy at the 
site strain) and slip interaction (given by the chemical misfit 
multiplied by the slip energy of the site) 
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This misfit approximation requires significantly less computing 
resources to determine compared with direct calculations. 
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Figure 3 : Mg basal screw dislocation core geometry with atomically resolved 
slip energy (top) and volumetric strain (bottom). The equilibrium geometry is 
found using first-principles flexible boundary condition methods. The circles 
show positions of atomic rows repeated out of the page (the [21 10] direction), 
while the magnitude of arrows between neighboring sites represent the relative 
displacement of the neighboring rows out of the page. The arrow magnitude is 
maximum for the partial slip. The dislocation splits into two partials — whose 
cores are centered on a closed circuit of arrows — with screw components (dis- 
placements out of the page) of 0.5a, and edge components (displacements in the 
horizontal slip plane) of 0.289a, separated by 3.9a. The slip energy is computed 
from the atomic geometry based on the relative displacements of neighbors, 
and has a maximum value of 20meV at the green sites. The volumetric strain 
is computed from the average change in nearest neighbor distances, and has a 
maximum tensile value of +4.6% at the cyan sites and maximum compressive 
value of -2.7% at the magenta sites. 



3. Dislocation / solute interactions 

Figure [3] and Figure [4] show the first-principles equilibrium 
basal screw and edge dislocation core geometries analyzed in 
terms of size and slip at each atomic site. The screw disloca- 
tion has displacements (slip) parallel to the line direction (out 
of the page), and the edge dislocation has slip perpendicular to 
the line direction, representing the two limiting cases for basal 
dislocation geometry. The equilibrium screw dislocation geom- 
etry dissociates into |[1010] and | [01 10] partial dislocations 
separated by 3.9a of 12 stacking fault, where a is the basal lat- 
tice spacing. The equilibrium edge dislocation also dissociates 
into partial dislocations (half the slip of the full dislocation), 
separated by 6.7a of stacking fault. The "cores" of the partials 
have the largest change in local geometry from bond-stretching 
and bond-bending. For a screw dislocation, there is no far-field 
dilation (tension or compression); however, we find a large vol- 
umetric strain in the partial cores. The strain is comparable to 
that in the edge dislocation core, which does have a far-field 
dilation strain. We quantify the bond-bending with the relative 
displacement d of each atomic row, and average the energy of 
stacking faults with displacement d of the six neighbors to give 
a slip energy E s n p . This energy is localized to the partial cores 
and the stacking fault separating them. A solute's binding en- 
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Figure 4: Mg basal edge dislocation core geometry with atomically resolved 
slip energy (top) and volumetric strain (bottom). The equilibrium geometry is 
found using first-principles flexible boundary condition methods. The circles 
show positions of atomic rows repeated out of the page (the [01 10] direction), 
while the magnitude of arrows between neighboring sites represent the relative 
displacement of the neighboring rows along the horizontal slip plane. The ar- 
row magnitude is maximum for the partial slip. The dislocation splits into two 
partials — whose cores are centered on a closed circuit of arrows — with edge 
components (displacements in the horizontal slip plane) of 0.5a, and screw 
components (displacements out of the page) of 0.289a, separated by 6.7a. The 
slip energy is computed from the atomic geometry based on the relative dis- 
placements of neighbors, and has a maximum value of 25meV at the green 
sites. The volumetric strain is computed from the average change in nearest 
neighbor distances, and has a maximum tensile value of +5.3% at the cyan 
sites and maximum compressive value of -4.9% at the magenta sites. 



ergy will change in the dislocations due to the different local 
geometry; moreover, as the cores show the greatest change in 
local environment, we expect solutes to provide the strongest 
interactions from sites within the dislocation cores. 

Figure B] shows the first-principles interaction energy — the 
change in the solute binding energy — for an Al solute with 
screw and edge dislocation cores. We find the strongest inter- 
action at the site of maximum compression near the center of 
the dislocation partial cores for both the edge and the screw ge- 
ometries. We expect this, since Al is smaller than Mg and is 
affected by the size changes in the core. The interaction force 
(the derivative of the interaction energy in the slip direction) is 
largest between the partial cores and the stacking fault for both 
dislocations, with nearly equal magnitudes for edge and screw. 
The similar values for the screw and edge forces are surprising 
as elasticity theory predicts a very weak far-field interaction 
for the screw dislocation! 12 1. If we a Pply anisotropic elastic- 
ity theory at a distance of c/4 from the slip plane (where c/2 
is the distance between parallel basal slip planes) to the screw 
and edge dislocations separated into mixed partials (|[1100] 
and | [1010]), anisotropic elasticity theory predicts maximum 
and minimum volumetric strains in the partial cores of ±7.2% 
(screw) and ±12.7% (edge). This is significantly different from 
our atomic-scale values of -2.7%, +4.6% (screw) and -4.9%, 
+5.3% (edge). Moreover, elasticity theory predicts the maxi- 
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Figure 5 : Al interaction energies with Mg screw (top) and edge (bottom) dis- 
locations, from first principles. We substitute Al atoms at different sites in and 
around the partial cores, and compute the energy differences after relaxation. 
The cyan sites show points of maximum binding energy while the magenta sites 
show maximum interaction force — i.e., greatest change in binding energy along 
the slip plane (c.f. Tablefllfor numerical values). The sites between the partial 
cores and stacking faults have the largest interaction force for both dislocation 
types. The strong Al interaction with the screw dislocation is surprising, and 
can only be accurately resolved through the atomic-scale calculation of the full 
dislocation core geometry. 



mum interaction forces for solutes to be the maximum interac- 
tion energy divided by 2A, and that the solute/screw interaction 
force will be 58% of the solute/edge interaction force. All of the 
elasticity predictions contradict our atomic-scale calculations: 
equal solute/screw and solute/edge interaction forces, energy to 
force ratios of 6-8 A, and smaller and unbalanced maximum and 
minimum volumetric strains. All of these differences highlight 
the sensitivity of the solute-dislocation interaction to the partial 
core geometries and the need for a first-principles approach to 
compute stress-free atomic-scale dislocation core geometries. 

We replace the elasticity approach to solute-dislocation in- 
teraction with a fully atomic-scale approach based on size and 
chemical misfits. The long-range strain and stress fields of a 
dislocation are known from anisotropic elasticity [2 1|, and so 
the change in local volume around a solute (size misfit Sb, the 
logarithmic derivative of the lattice constant with solute con- 
centration) controls the solute-dislocation interaction. In the 
dislocation core, the interaction energy is largest, but it is no 
longer described by elasticity. Despite this, we can use the size 
misfit even in the partial cores by computing the change in local 
volume for each atomic site. This approximation is the largest 
contribution to the interaction energy; the next largest contribu- 
tion is from the slip in the partial cores and the stacking fault. 
In the same way that we determine a slip energy in the core 
of the dislocation, we find that solutes change the response of 
the crystal to slip (chemical misfit ssfe, the logarithmic deriva- 



tive of the stacking-fault energy with solute concentration). The 
chemical misfit determines how the atomically-resolved slip en- 
ergy (c.f., Figure[3]and FigureRl will change with the addition 
of a solute. Adding this contribution to the size misfit produces 
an accurate approximation of the solute-dislocation interaction 
energy using the dislocation core geometry. 



Table I : Al interaction energy with Mg basal screw and edge dislocations from 
direct solute substitution calculations (c.f. Figure]?) and misfit approximations 
based on the dislocation geometry. The maximum binding energy (meV) and 
maximum interaction force (meV/A) determine the attraction of solutes to dis- 
locations and solid-solution strengthening. The misfit approximation uses two 
misfits with the atomic-scale dislocation geometry: change in local volume, 
and bond-bending from slip to compute interaction energies for solute with 
the dislocation core. For both dislocation geometries, the misfit approximation 
correctly captures the interaction energies and forces compared with the more 
computationally intensive direct calculation. 
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Table [T] compares the solute-dislocation interaction for Al 
computed using direct substitution into the dislocation cores 
with our combined size and chemical misfit approach. The re- 
sponse of Al to changes in local volume and slip capture most of 
the interaction energy and forces in the dislocation core. More- 
over, the strain and slip energies are taken directly from the 
equilibrium core geometries, and the size and chemical misfit 
calculations include the local response of Mg atoms neighbor- 
ing the Al solute. The size-misfit approximation is expected to 
be accurate in the far-field, and is accurate even in the partial 
cores where the interaction is the strongest. The slip energy is 
needed to represent the stacking fault region between the par- 
tials, and the center of the partial cores themselves. The maxi- 
mum interaction forces from the misfit calculation are accurate 
to within 5% of the direct calculations for both screw and edge 
geometries, with deviations of only 5meV for the interaction 
energies; hence, we can predict the interaction energies of other 
solutes by using the size and chemical misfits combined with 
the equilibrium dislocation core geometries. This permits us to 
use much simpler and computationally efficient first-principles 
calculations of misfits with our first-principles calculation of 
pure Mg dislocation cores. 

Figure [6] provides the formulae for the maximum solute- 
dislocation interaction force in terms of misfits. The interaction 
force from the first-principles solute binding energy is F max = 
max |m-VZ?bindingl where m is the dislocation slip direction ([1 100] 
for screw and [1 120] for edge), and ^binding is the sum of the so- 
lute slip and size interaction energies from EquationlT] The gra- 
dients of the solute size and slip interactions are calculated as a 
centered difference in the slip direction for sites near the partial 
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Figure 6: Formulae for computation of maximum solute-dislocation interaction 
forces from £;, and ssfe for edge (top) and screw (bottom) dislocations. Each 
solute substitutional site in the screw or edge dislocation has a local volumetric 
strain and slip energy; we use the size and chemical misfits to approximate the 
solute-dislocation interaction energy by linearly scaling those energies. The 
derivative of the solute-dislocation interaction energy along the slip direction 
gives the interaction force, and the maximum interaction force appears as input 
in our strengthening model (Equation |10) . The particular site that produces the 
largest interaction force changes depending on the magnitudes of the misfits £;, 
and £sfe, leading to different regimes. The slopes of the lines separating the 
different regimes are indicated in red. 



dislocation cores where the gradient is the largest. The gradi- 
ents were scaled by the size and chemical misfits for each so- 
lute and summed to determine the maximum interaction force. 
The end result is that, for any combination of ej, and esfe, the 
strongest interaction force site is known, and its pinning force 
can be computed as a linear combination of the two parameter; 
however, the coefficients change as different magnitudes of £/, 
and £sfe will select different sites as the strongest interaction 
force site. Figure [6] graphically summarizes all of the formulae 
for each range of s/, and ssfe- For example, Al has Sb - -0. 1 15 
and ssfe = -1-252; then, for the edge dislocation, the inter- 
action force is given by the formula in the southwest corner 
171.3s/, + 2.62ssfeI = 11.5meV/A, and for the screw disloca- 
tion, the interaction force is given by the formula in the nar- 
row southwest wedge |102e/, - 0.370ssfeI = 11.2meV/A. For 
Zn, ei, - -0.153 and ssfe = +0.317, so the edge dislocation 
interaction force is given by the formula in the south section 
|101efe - 0.412£sfeI = 15.6meV/A and for the screw disloca- 
tion, the interaction force is given by the formula in the south 
section |104e fe - 0.599s S feI = 16.1meV/A. 

4. Solid-solution strengthening model 

To predict solid-solution strengthening, we use our first- 
principles atomic-scale solute-dislocation interaction calcula- 



tion as input to a dilute-concentration, weak-obstacle model for 
solid-solution strengthening from Fleischer|23|. As the disloc- 
ation moves in the slip plane under stress, it encounters ran- 
domly placed immobile solute atoms each of which provides a 
"pinning" force up to the maximum solute-dislocation interac- 
tion force F max . This point-pinning model is applicable for iso- 
lated obstacles with a short-ranged — on the scale of the solute- 
separation distance — interaction force between solute and dis- 
locations, and hence will capture the dilute-concentration limit. 
It also assumes a random solute distribution that does not re- 
arrange due to the dislocation strain-fields; hence, it is opera- 
ble at temperatures where no appreciable solute diffusion oc- 
curs. The nearly straight dislocation bows at the solute un- 
til it reaches a critical bowing angle C and the line tension <5 
pulls the dislocation past the solute. At the critical angle 9 C , 
Fmax - 2£sin((7r - f? c )/2). For weak obstacles, the dislocation 
is nearly straight, and the critical bowing angle is only slightly 
smaller than n, hence n - 6 C a F max /£ <s 1. The mean dis- 
tance between pinning points, L, is the average distance be- 
tween nearest randomly placed solutes in a circular wedge of 
angle a — (n - G c )/2, where L - Jn/(2p s a), and p s is the den- 
sity of solutes in the slip plane. The solute density is 2c s (atomic 
concentration c s ) per y/3b 2 /2 (Burgers vector b — a, basal lat- 
tice constant), assuming that the solute can appear either above 
or below the slip plane; hence, 



Hence we need to consider the contributions to strength 
from both edge and screw dislocations for all solutes in our 
modified-Fleischer model. The concentration-independent pref- 
actor in Equation 10 is the solute "strengthening potency." The 
solubility of a solute controls the maximum possible yjcl, and 
hence the maximum possible strengthening of a solute in an 
alloy. Combining the calculation of interaction strengths from 
their misfits, the strength potency can be fit to an approximate 
functional form 



Atcrsscoood ~ (38.9MPa) 



\\0.176/ 



1 5.67/ 



•,3/2 
g A g SFE I r/2 

2.98 J ' Cs 

(11) 



The negative £/,£sfe term shows that the effect of the size and 
chemical misfits is not completely independent, but that due to 
the geometry of the dislocation core — where the sites of largest 
slip energy are neighboring sites of largest strain energy — a 
larger gradient can be produced by having the size and slip mis- 
fits with opposite signs. However, the magnitude of the denom- 
inators also shows that the size misfit is the dominant interac- 
tion. The fit had a maximum error of 10%, with the average 
error being 3%. As this model ignores thermally activated pro- 
cesses to overcome solute obstacles, it is appropriate at temper- 
atures without appreciable diffusion, but the fundamental inter- 
action data is usable for high temperature strengthening models 



1/2 = ( n^b 2 V /2 = (3n 2 y' 4 b I 6 \ 1/2 

P,(*-6 e )) ~U(7T-0 C ) C J 2 Umax/ ^ 

(8) 
That is, the distance between pinning points increases with de- 
creasing solute concentration c s and with decreasing interac- 
tion force. Given the mean solute spacing, the strengthening 
(increase in critical-resolved shear stress in the basal plane nec- 
essary to overcome the solute restraining force) is: 



as well. All of the material information that enters Equation 10 



-1/2 



comes completely from first principles: crystal structure, lattice 
and elastic constants, dislocation core geometries, and solute 
misfits. 
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which scales as F m ' ax and -\/c^. This is the change for a (nearly) 
straight dislocation line of a given character ranging from edge 
to screw; the interaction energy and line tension correspond to 
the particular line orientation. The strengthening of an aver- 
age dislocation loop — which continuously ranges from edge to 
screw orientations — is weighted by the line tension. The shape 
of a basal loop can be estimated as an ellipse with axial ratio 
a = K e /K s , where K e = 25.6GPa and K s = 18.6GPa are the 
edge and screw line energy prefactors[20| as determined from 
the first-principles elastic constants Cy. For Mg, this gives a 
weight of 0.539 for the screw dislocation and 0.461 for the 
edge. The screw and edge line tensions are S s - \K s b 2 - 
591meV/A and 8 e = \K e b 2 = 817meV/A; so, the dislocation 
loop weighted-average strengthening for each solute is 
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where AtcrsS(0001) is i n MPa and the interaction forces are in 

meV/A. 
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Figure 7: Solid-solution strengthening potency (contours) versus their size (ej,) 
and chemical (esfe) misfits for 29 different solutes in Mg. The misfits, com- 
bined with the screw and edge dislocation geometries (c.f. FigurefSland Fig- 
urepl determine the maximum interaction forces. In the dilute limit for weak 
obstacles, the strengthening — characterized by a change in the basal critical- 
resolved shear stress Atcrss(OOOI) — scales with ^[cl and with the interaction 
forces to the 3/2 power in Equation |10| We can efficiently predict the solute 
strengthening for a whole range of solutes from the first-principles core geom- 
etry. This "design map" shows solutes with comparable strengthening potency, 
and also gives a guide to solubility — which decreases with increasing misfit 
magnitudes. 




FigurelTlcombines the first-principles data for size and chem- 
ical misfits with the simple solid-solution strengthening model 



of Equation 10 into a "design map" for the strengthening poten- 



cies of 29 different solutes. The contours show equal strength- 
ening potency versus size and chemical misfits. The size and 
chemical misfits can be easily computed for any substitutional 
solute in the periodic table, and then the screw and edge dis- 
location maximum interaction forces to predict the change in 



critical-resolved shear stress for basal slip from Equation 10 



This map provides a rational method to select equipotent solutes 
to replace less "desirable" elements based on high mass, cost, 
or other processing concerns. In addition, the misfits also give 
rough qualitative information about solubility, as larger mis- 
fit magnitudes lead to lower solubilities. For example, for Zn 
with a potency Atcrss(OOOI)/ sfcs - 32.5MPa, the correspond- 
ing small magnitudes of £/, and ssfe suggest a high solubility. 
Alternatively, the much larger potencies of Ir (172MPa) and K 
(161MPa) are due to the larger magnitudes of £/, and ssfe, giv- 
ing rise to a much lower solubility. Yttrium represents a rea- 
sonable compromise between strength and solubility. 

Table [2] provides misfits and potencies for all 29 solutes, 
as well as comparisons to available experimental data (for Al, 
Zn, Bi, Cd, In, Li, Pb, Sn, and Tl). For comparison, we look 
to single-crystal, low temperature measurements of critical re- 
solved shear stress of Mg in the low concentration limit. For 
example, the available experimental data for potencies of com- 
mon solutes Al and Zn in Mg extrapolated to 0K[24, 29 1 give 
21.2MPa and 31MPa, versus our first-principles prediction of 
19.5MPa and 32.5MPa — validating our computation and mod- 
eling approach. For Bi and Pb, only plateau stress data is available — 
this is the strengthening divided by y/c^ at larger concentrations 
than considered here. We expect the plateau strengthening co- 
efficient to be a lower limit on our dilute concentration potency, 
and so our data remains consistent with the available experi- 
mental data. 

5. Conclusions 

The basal-slip strengthening design map for Mg represents 
an important development in predicting chemical effects on strength 
ening in an accurate and efficient manner. In addition to pre- 
dicting strength, we find a strong solute interaction with screw 
dislocations equal to that of edge dislocations: this demon- 
strates the need for first-principles quantum-mechanical calc- 
ulations with flexible boundary conditions to reveal defect in- 
teractions. Maximizing strength requires careful consideration 
of the important tradeoff between high strengthening capac- 
ity (i.e., large solute misfits) and high solubility (low misfit 
magnitudes). Predicting strength above dilute concentrations, 
at higher temperatures and including other mechanisms such 
as forest strengthening! 30 1 will require the fundamental input 
from the present solute strengthening predictions. Moreover, 
cross-slip in Mg depends on the constriction of basal screw 
dislocations, and should be aided by solutes with positive ssfe 
to increase the basal stacking-fault energy. By finding solutes 
that can improve cross-slip while simultaneously strengthen- 
ing basal slip — as approached here — the strength anisotropy of 



Mg alloys can be reduced, improving ductility. The natural 
extension of the computational methodology developed herein 
is to predict solute effects on slip in prismatic planes, includ- 
ing increasing cross-slip for ductility enhancement. In addi- 
tion, atomically-resolved strengthening predictions hold great 
promise for the design of other technologically important HCP 
metals such as Ti and Zr. 
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